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Abstract 

We point out under the assumption of Majorana neutrinos that a lower bound 
on the MNS matrix element |^7e3p can be derived by using constraint imposed 
by neutrinoless double beta decay experiments and by positive detection of 
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experiments, {m)pp/ {m)p. It means that once the neutrino mass is detected 
in the bound-sensitive region one must soon observe signal in neutrinoless 
double beta decay experiments. 
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I. INTRODUCTION 



There exist, by now, accumulated evidences in the atmospheric and the solar neutrino 
observations that neutrinos do oscillate. The existence of neutrino oscillation is further 
strengthened by the result of the first long-baseline man-made beam experiment K2K, in 
particular by their latest result . They constitute the first compelling evidence for physics 
beyond the standard model of particle physics. In particular, we have leaned that an almost 
maximal mixing angle ^23 is required to account for the atmospheric neutrino anomaly, which 
is quite unexpected from our experience in the quark sector. 

What is even more surprizing to us is that, according to the latest global analyses of the 
current solar neutrino data the angle 612 which is responsible for the solar neutrino 

oscillation is likely to be large though may not necessarily be maximal. It is in sharp contrast 
with the fact that the remaining mixing angle 613 is constrained to be small, s^3<0.03, by 
the reactor experiments [0]. Given the current status of our understanding of the structure 



of lepton flavor mixing matrix, the MNS matrix [jTO| , it would be nice if there are any hints 
on how small is the angle 6'i3.|^ 

In this paper, we try to pursue such a possibility and point out that one can derive a 
lower bound on sf^ = |f/e3p through joint efforts by double beta decay experiments and by 
direct mass determination either by single beta decay or by cosmological observations. We 
assume in this paper that neutrinos are Majorana particle to rely on the bound imposed by 
double beta decay experiments. 

Before getting into the bussiness, let us briefly summarize existing knowledge on how 
^13 can be measured, or further constrained. Most optimistically, the next generation long 
baseline experiments, MINOS ||12[, JHF [|13|, and OPERA |]14| will observe z/e appearance 
events and measure the angle 613. Most notably, the JHF can probe sin^26'i3> several 



xlO~^ in its phase I ||T3|. If realized, a large- volume reactor experiment [15] can also probe 
the similar (to a slightly shallower) region. If the angle is smaller than the sensitivity region 
of these experiments, we have to wait for future supermassive detector experiments, utilizing 



either low energy conventional superbeams or neutrino factoties. (See e.g., for references 
cited therein.) If nature is so unkind as to tune the angle extremely small, sin^ 26'i3 10~^, 



then the only way to detect its effect would be via supernova neutrinos [17 



See e.g., [11 1 for a summary of remaining issues in three flavor mixing scheme of neutrinos. 
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II. CONSTRAINT FROM NEUTRINOLESS DOUBLE BETA DECAY 



Let us start by examining constraint from double beta decay. We use throughout this 
paper the standard notation of the MNS matrix: 



U 



C12C13 



■S12C13 



~'5l2C23 — Ci2S23Si3e*'^ C12C23 — 5125235136**^ 523C13 
512^23 ~ Ci2C23Si3e*'^ " C12S23 — Si2C23Si3e*'' C23C13 



(1) 



Using the notation, the observable in neutrinoless double beta decay experiments can be 
expressed as 

3 



(m) 



i=l 



(2) 



where (i=l, 2, 3) denote neutrino mass eigenvalues, Uei are the elements in the first 
low of the MNS matrix, and /? and 7 are the extra CP-violating phases characteristic to 
Majorana neutrinos ||18|,|19[. We have used in the second line of (S) the Majorana phases in 



the convention of |20|. There have been large number of papers quite recently which devoted 
to extract constraints from neutrinoless double beta decay experiments p0|-p2|. 



We define the neutrino mass-squared difference as Am?- 



m\ in this paper. In 



the following analysis, we must distinguish the two different neutrino mass patterns, the 
normal (Am23 > 0) vs. inverted (Am23 < 0) mass hierarchies. We use the convention that 
m3 is the largest (smallest) mass in the normal (inverted) mass hierarchy so that the angles 
Q\2 and 6'23 are always responsible for the solar and the atmospheric neutrino oscillations, 
respectively. We therefore sometimes use the notations Am23 = Am^^^ and Am^2 = Amg 
to emphasize that they are experimentally (the latter to be) measured quantities. Because 
of the hierarchy of mass scales, Amg/Am^j^ ^ 1, Am^2 can be made always positive as 
far as ^12 is taken in its full range [0, 71/2] ||23|| . 

In order to derive constraint on mixing parameters we need the classification. 
Case A: 



miCi2Ci3e "-^ 77123^20136+*^ 



> 



m3Si3 



(3) 



Case B: 



miCi2Ci3e *^ -I- m2Si2Ci3e"'"*^ 
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(4) 
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Both types of mass hierarchies are allowed in the cases A and B. For a given experimental 
upper bound on {m)pj3, one can derive a lower (upper) bound on 5^3 in the case A (B). 
Therefore, we start with the case A. 



A. Case A 



In this case, the lower bound on {m)/3/3 can be obtained as in the following way; 



{m)fifs > (miCi2 + msSia) cos/3 - i{mic^2 - "^2^12) sin/3 



= c\.^\J rn\c\2 + '^2^12 + 2mim2cf2Si2 cos 2/3 — m2,s\^. (5) 

Noticing that the right-hand-side of (^ has a minimum at cos 2/3 = —1, we obtain the 
inequality 

{m)pp > miCi2 - ^2^12 - m3Si3. (6) 

If a neutrinoless double beta decay experiment imposes the bound {m) < {rnYpp the 
lower bound on 5^3 results; 



\mlc\2 — 1TI2S12 




/ \ exp 




micf2 - ^25^2 





4s > ' ,1 .2 1 (7) 



We rewrite the lower bound into the one expressed by the heaviest neutrino mass mn 
and Am^ measured by the atmospheric and the solar neutrino experiments. We note that 
the neutrino masses rrii (i=l,2,3) can be parametrized by Am^^^ = Am^s, AitLq = Amf2, 
and a absolute mass scale tuh- We take mu = and itih = "^2 for the normal and the 
inverted mass hierarchies, respectively. Then, 



mi = y mjj - Amlt^ - Ami, ^2 = y mj^ - Amlt^, mg = m^, (8) 
for the normal mass hierarchy, and 



rui = \J rn?H - Am%, m2 = mH, nis = sjmjj - \Amlt^\, (9) 

for the inverted mass hierarchy. 

We express the lower bound on 5^3 by using the dimensionless ratios 



rriH ttlh ttlh 

The lower bound on s\.^ reads for each type of mass hierarchy as follows: 
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(i) Normal mass hierarchy; mn = fn^ 



4> 



2 ^2 



atm 



atm'-' 12 



R 



P(3 



atm 



^©^12 



atm.^12 



+ 1 



(11) 



(ii) Inverted mass hierarchy; itlh = Tn,2 



4> 





\/l-4c?2- 


s2 
■'12 






/i^2 --i2 f^2 
'©'-12 "'12 


+ 





(12) 



The parameter is extremely small, Tq = ArriQ/m'jj < 10~^ for the possible best sensi- 
tivity of ~ 0.3 eV. (See later.) Hence, it is an excellent approximation to ignore Tq unless 
the sensitivity goes down to very close to Am; 



2 

atm 



SO that it is comparable with 1 — r^^^ in 



the normal mass hierarchy case. Ignoring the lower bound on sfn greatly simplifies: 



4> 



cos 26^12 1 



1 — ^Itm — Rpli 



4> 



COS 26*12 1 
I cos 20 



/V.2 
' atm 



121 



R 



cos 2^12 I + \l 



(normal mass hierarchy), 
(inverted mass hierarchy). 



(13) 
(14) 



atm 



We finally note that in the degenerate mass limit Tatm —>■ the lower bound ceases to 
distinguish between the mass patterns, and has a universal form 

(degenerate mass limit). (15) 



COS 26^12 




- Rf3(3 




cos 26*12 


+ 1 



B. Case B 



For completeness, we treat the case B, which yields the upper bound on sfg. Proceeding 



via the similar way toward (^) we obtain the inequality 



-13 



miCi2 + rra2Si2 



(16) 



which is saturated at cos 2/5 = +1. Then, the upper bound on s^g for a given experimental 
bound on (m)^^ entails after ignoring as 



43< 



4< 



atm 



+ R 



Pf3 



atm 



+ 1 



l + R 



(normal mass hierarchy). 



(inverted mass hierarchy). 



(17) 



(18) 



atm 



It can give a stronger bound than the CHOOZ limit only for the normal hierarchy and if 
Rjsp < 0.03. We do not discuss it further because it is outside of our analysis in section IV. 



5 



III. CONNECTION BETWEEN Mh AND THE OBSERVABLE IN DIRECT MASS 

MEASUREMENTS 



Before entering into the actual analysis of the bound, we clarify the relationship between 
the largest mass mn (= "^s for normal, and = m2 for inverted mass hierarchies) and the 
observable in direct mass measurement in single beta decay experiments.^] We argue that itih 
can be identified as the observable in such experiments in a good approximation. Suppose 
that the neutrino masses are hierarchical and obey either ms ^ 1112 — nii (normal mass 
hierarchy), or m2 ^ mi ^ m-^ (inverted mass hierarchy). In the former case, it is obvious 



that rriz = mn is the observable. In the latter case, it was shown in [^] that the observable 
in direct mass measurements (m)^ is given by 



= I.. 10 (19) 







2 






2 



where n is the dimension of the subspace of (approximately) degenerate mass neutrinos, 
and 72 = 2 in the case under discussion. Then, (m)/3 = m2 = mn in a good approximation. 
In the opposite extreme, mf ^ Am^j^ which is usually referred to as "almost degenerate 



neutrinos", (|T9D with n = 3 tells us that (m)^ = m^^ in a very good approximation. 

Thus, (m)^ = tuh holds in both extreme, hierarchical and degenerate mass neutrinos. 
This discussion strongly suggests that (m) ^ is reasonably well approximated by itih even in 
the intermediate region. 



IV. ANALYSIS OF THE LOWER BOUND ON Sfg 

We present in Fig.l the lower bound on 5^3 as a function of mu for various values of 
cos2^^i2; the region lower-right to each curve is excluded. The shaded area indicates the 
region excluded by the CHOOZ experiment 0, which we approximate as 5^3 < 0.03.0 We 
take two typical values of (m)^, 0.34 eV and 0.1 eV. We present in Fig. 1 only the results 



For the present status and the future prospects of this type of experiments, see e.g., the website 



of a recent conference devoted to the topics |24|. 

^ While the precise value of the CHOOZ constraint actually depends upon the value of |Am^j^| 
P], we do not elaborate this point in this paper. 
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for the normal mass hierarchy. It is because the results barely changes in the case of inverted 
mass hierarchy; there is practically no diference in (a) (^)^^^ = 0.34 eV case, and the curve 
shifts toward the right by about 4 % in (b) {fny^j^ = 0.1 eV case at the best fit value of 
LMA solution (see below). 

We note that the former value of {mYp^ corresponds to the present 90 % CL bound by 
Heidelberg- Moscow group p6[ , while the latter indicates a modest sensitivity to be achieved 



in near future by CUORE GENIUS |28|, and by MOON [||] double beta decay exper- 
iments. We would like to remind the readers that the present strongest upper limit on (m)^ 



is from the Mainz collaboration |30|, (m)^ < 2.2 eV (95 % CL). A similar bound (m)^ < 2.5 
eV (95 % CL) is derived by the Troitsk group [^. The sensitivity of the proposed KATRIN 



experiment is expected to extend to {'m)j3 < 0.3 eV p2 



In the numerical analysis in this section we assume \Aml^J = 3 X 10~^ eV^, the best fit 
value ||3^ of the combined data of Super-Kamiokande [|l| and the K2K experiments. We 
ignore Am^ in most part of our analysis and its effect is detectable only in a limited region 
outside of the sensitivity of the KATRIN experiment, as shown in Fig. 2. 

We note that the allowed region of the mixing angle 612 has a large uncertainty. In 
particular, its value at the largest end is far more uncertain compared with the accuracy 
we need to determine where is the bound-sensitive region. We just quote for illustration 
the allowed region for the LMA solution given in 0; ~ 0.67 > cos 26*12 > 0.19 (95 % CL), 
> 0.099 (99 % CL), and > -0.024 (99.73 % CL). Then, the allowed region of the LMA 
solution extends even at 95 % CL to the right-most contour in Fig. 1. Therefore, the Mainz 
and the Troitsk experiments begin to touch the parameter region already at their present 
sensitivities. Similarly, the 95 % allowed regions of the LOW and the VAC solutions extend, 
very roughly to, 0.30 > cos 2^12 > 0.026 and 0.57 > | cos 2^12! > 0.30, respectively 0. 

We observe in Fig. lb that for {mY^p = 0.1 eV the KATRIN experiment start to 
cover a large portion of the parameter region of the LMA solution toward its best fit value 
of cos 26*12 = 0.48 (0.46 in 0). Therefore, if the LMA solution is confirmed by the 
KamLAND reactor experiment and if the KATRIN experiment detects direct neutrino mass, 
we would hope that we will have a lower bound on sf^ at the same time. 

In Fig. 2 we present the similar plot as in Fig. 1 but for {m)'^^^ = 0.01 eV. In this case 
the cases of the normal vs. the inverted mass hierarchies differ clearly from each other, as 



one can see by comparing Figs. 2a and 2b. In region of mn ~ y |Am^j„| ~ 0.06 eV the lines 
of positive and negative cos 26*12 start to split by nonvanishing correction, as indicated 
by the dashed lines in Fig. 2. 
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The significant feature of the lower bound we have obtained is that it exists in a narrow 
window of mn, the highest neutrino mass, whose values depend very sensitively on {m)pp, 
the observable in neutrinoless double beta decay experiment. We now try to characterize 
the region, the bound-sensitive region, as a function of the experimental upper limit {jnYp^ 
and {m)p. In the degenerate mass limit Tatm 0, it can readily be done by using (pISl). By 
demanding the consistency with the CHOOZ bound, we obtain 

0.971 cos 2^12 1 - 0.03 < Rag ^ \ < I cos2^i2|. (20) 

Therefore, the bound-sensitive region is characterized by the ratio of these two experimental 
observables in the degenerate mass approximation. This feature prevails to certain extent 
beyond the approximation as we will see. 

In Fig. 3 we present plots of the bound- sensitive region expressed by the ratio Rf^jj = 
{mYp^/{m)j3 as a function of cos26'i2 for the normal and the inverted mass hierarchies in 
regions for (a) LMA and (b) LOW solar neutrino solutions. The region covered is determined 
as the best fit value ± 10 % in tan^ 9i2. For the LMA solutions it roughly corresponds to the 
accuracy to be achieved in KamLAND ||3^. We do not present the case of vacuum solution 
because the relevant region roughly overlaps with that of LMA solution. In the degenerate 
mass limit we must have a universal curve (|20|) and the splitting between the normal and 
the inverted hierarchy cases in Fig. 3 represents corrections by the effect of nonvanishing 
Tatm- We do not plot the scaling curve ( pOD because it is virtually identical with the one for 
the inverted mass hierarchy. 

Notice that the ratio Rpp can be small for each mass pattern because of the suppression of 
contribution from highest mass neutrinos by small in the case of normal mass hierarchy, 
and by possible cancellation in the inverted mass hierarchy. Once Rp^^ is given, it is easy for 
the readers to read off the lower bound on 5^3 for a given value of 9i2 because the function of 
lower bound (p!3D and ([T^) are approximately linear within the narrow strip < sfg < 0.03. 

We note, however, that several highly nontrivial requirements must be met in order 
to extract the lowest value of 5^3. We see from fig. 3 that the accuracy of direct mass 
measurement of {m)p at given {mY^I^ must be better than roughly ± 5 % for the LMA 
and ± 20 % for the LOW solutions, respectively. Furthermore, the precise values of mixing 
parameters, in particular, cos 2^12 must be known to better than ± 10 % level. 

Finally, some remarks are in order: 

(1) Suppose that in a future time the solar mixing angle 612 is determined by some ingenious 
experiments. Assume then that we have obtained the lower bound on by KATRIN 
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observation of neutrino mass at the time of sensitivity {mY^p =0.1 eV of double beta decay 
experiment. Now, further suppose that the latter experiment would have been improved so 
that {mYpp <^ 0.1 eV. Figure 3 tells us that then we are in the excluded region. What does 
it mean? 

The answer is; it means that the double beta decay experiment must see positive events 
before going down to (m)^^^ <^ 0.1 eV. If not, it is the indication that nature chooses a 
different lepton flavor mixing scheme from the standard one, or our assumption of Majorana 
neutrinos is wrong. The last possibility is a rare case of excluding Majorana hypothesis in 
spite of finite resolution of {m)j3p measurement, which takes place owing to the CHOOZ 
constraint. It is very exciting that a discovery in an experiment either signals discovery 
of another experiment, or indicate radically different views of lepton sector from the one 
useally accepted. 

(2) We have explored in this paper features of the lower bound on sf^ which can be derived 
by using constraints imposed by the single and the double beta decay experiments. The 
similar consideration can be done to constrain different parameters such as solar mixing 



angle 9i2, for example. A more generic analysis of the bound is in progress |^ 



(3) While we have focused on beta decay experiments in this paper, the possibility of obtain- 
ing better knowledge of neutrino mass in terms of cosmological observation must be pursuit. 
See e.g., [pB . 
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FIG. 1. The lower bound on 5^3 is displayed as a function of mn for various values of | cos 2Q\2\ 
as indicated in the figure for (a) {mf^^ = 0.34 eV and (b) = 0.1 eV. The region lower-right to 

each curve is excluded. The shaded area indicates the region excluded by the CHOOZ experiment. 
Only the case of normal mass hierarchy is presented. (See the text.) 
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FIG. 2. The same as in Fig. 1 but with {m)^^^ = 0.01 eV. Figures 2a and 2b are for cases 
of the normal and the inverted mass hierarchies, respectively. The upper and lower dashed lines 
around curves of | cos 2^12! = 0.2 and 0.15 in both figures are for cos 2^12 < and cos 2^12 > 
cases, respectively, with Auiq = 4.8 x 10~^ eV. 
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FIG. 3. The region of i?^^ = {fnY^^ / {m) p in which the lower bound exists is exhibited as 
a region between two Hnes as a function of | cos 2^i2 1 . The sohd and dotted hues correspond to 
the normal and the inverted mass hierarchies, respectively. Figures 3a and 3b cover roughly the 
parameter regions of the LMA and the LOW solar neutrino solutions, respectively. 
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